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^ \ Abstract. Gravitational microlensing is a new technique for studying 

the surfaces of distant stars. A point mass lens, usually a low-mass star 
\ from the disk, will typically resolve the surface of a red giant in the Galac- 

(T) ■ tic bulge, as well as amplify its brightness by a factor of 10 or more. Such 

CN [ events are now detected in real time and can be followed up with preci- 

sion photometry and spectroscopy. The motivation for stellar imaging via 
microlensing lies in its unique ability to provide center-to-limb variation 
■ of spectral diagnostic lines. Such variation maps into a variation with 

04 . depth in the stellar atmosphere — a unique and very valuable constraint 

CN ■ for modeling the structure of red giants' atmospheres. We illustrate the 

technique on a recent successful observational campaign and discuss the 
qq ' implications for current stellar models. 

On ■ 

Or 
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' Gravitational microlensing is a new technique for studying the surfaces of distant 

d . stars (Sasselov 1996, 1997a), which is observationally viable for objects in the 

Galactic bulge (Lennon et al. 1996, Alcock et al. 1997). Microlensing events 
•*h . have typical total durations of weeks to months. While most lenses are point 

K> ' masses, many (~ 20%) of the sources monitored in the Galactic bulge are red 

giants and supergiants — their angular radii are comparable to the Einstein radius 
of a sub-solar mass lens (in the range of 50/ias). In addition to being resolved 
by most lenses, the projected disks of bulge giants are large enough to make 
the probability for a lens transit very high — in fact, at least one such event 
(MACHO Alert 95-30) was very well observed recently (Alcock et al. 1997), and 
two near misses occured this summer (MACHO Alert 97-BLG-10 & 97-BLG-56). 

We study stars by the radiation emergent from their surface layers. While 
the theory of numerical radiative transfer has made great progress over the last 
decade, its application to stellar atmospheres has been, for the most part, limited 
to one-dimensional equilibrium models. One reason was, partly, because of tech- 
nology — we were unable to observe directly the surfaces of other stars. Stellar 
interferometry, eclipsing binaries, and Doppler imaging, provide valuable surface 
information, like limb darkening and stellar spots. However, by its nature, mi- 
crolensing imaging allows us to observe the center-to-limb variation in spectral 
lines (and their profiles). Such observations are invaluable for the construction 
of detailed models of the atmosphere, as it has been long known from the study 
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of the Sun. The solar disk brightness distribution as a function of wavelength 
has been mapped into the depth distribution of temperature and density, which 
allows a complete synthesis of the observed solar spectrum (Kurucz 1990). The 
situation with stars has been confined to very much lower accuracy; the data 
cannot be used to build a model atmosphere, as in the case of the Sun. 

The 2-D and 3-D structure of stellar atmospheres other than the Sun is not 
well known. Theoretical modelling relies heavily on the solar atmosphere model 
(Nordlund & Dravins 1990), but this situation is unfortunate for a number of 
reasons. There has been an increased demand for accurate stellar models from 
a number of fields, especially regarding red giants and supergiants. Our Sun is 
not a good zero point for most of these applications. 



2. Gravitational Microlensing as a Stellar Microscope 

The amplification of a point source by a point mass, M, depends only on their 
projected separation d (Paczyhski 1986), 

Atj\ 62 + 2 m 

^ = d(cP + 4)V* (D 

where d is expressed in units of the angular radius of the Einstein ring of the lens, 
Q E = ([4GM/c 2 ][Z) LS /Z) OLJ D os ]) 1 / 2 , and D QL ,D LS , and D os are the distances 
between the observer, lens and source. The total flux received from an extended 
source is therefore obtained by integration over its infinitesimal elements, 

F(t) = / S rdrB{r) / d9A(d), (2) 

JO JO 

where B(r) is the surface brightness profile of the source in the projected polar 
coordinates (r, 6) around its center. For simplicity, we assume that the stellar 
emission is circularly symmetric, although this is not a limitation (see Heyrovsky 
& Loeb 1997; Heyrovsky, Loeb, & Sasselov 1997), and denote the source radius 
by R s = R s x (6 E D os ). All projected length scales are normalized by the Einstein 
ring radius. Using the law of cosines 

d = \d - f | = (d§ - 2d r cos 9 + r 2 ) 1 ' 2 , (3) 

where do = (b 2 + v 2 t 2 ) x l 2 is the projected separation between the lens and the 
source center, b is the projected impact parameter of the source center, v is the 
transverse velocity of the lens in units of 9 E D QL , and t = is the midpoint time 
of the lensing event. For a circularly symmetric source, the resulting light curve 
is time-symmetric with F(t) = F(—t). 

As a result of this dependence on B(r), when the source is resolved, F(t) 
is wavelength dependent through B(r), despite the achromatic nature of grav- 
itational lensing. For example, in the optical continua, B{r) will reflect the 
wavelength dependence of limb darkening, hence broad-band photometry will 
show chromaticity due to that effect (Gould & Welch 1996; Heyrovsky, Loeb, & 
Sasselov 1997; Valls-Gabaud 1997). Narrow-band photometry could reveal res- 
onant line emission in cool giants in a much more dramatic ways, as shown by 
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Figure 1. Line profiles at three different viewing angles (character- 
ized by their cosines //) on the stellar disk from the non-LTE model of 
the M95-30 red giant (Sasselov 1997a). a) Three profiles of the hydro- 
gen Ha line from disk center at ^=0.98 (bottom), through [i=0.5, to 
very near the limb at n=0.02 (top), b) Same for the hydrogen H/3 line. 
Note the difference in profile shape between the two lines, as well as 
their change with jjl. 



Loeb h Sasselov (1995). And finally, high-resolution spectra will reveal spectral 
line profile changes with time, representing the subtle center-to-limb variation 
of different atomic and molecular transitions. 



3. Probing the Atmospheres of Red Giants 
3.1. Stellar Microscopy 

Stellar disks are projected hemispheres, which implies an axisymmetric variation 
of the intensity with position as it maps into a variation with depth. In general, 
the emergent continuum radiation at the center of the disk is formed deeper than 
the radiation we see near the limb. Over a very wide spectral range for stars of 
different temperature, the center-to-limb variation of the continuum emission is 
manifested as limb darkening. Unlike continuum radiation, spectral lines provide 
a larger choice of center-to-limb variations. In the majority of stars we deal with 
absorption lines of different strength, arising from atomic or molecular bound- 
bound transitions. Most moderately strong and weak lines diminish towards 
the limb (Figure 1). This variation could be observed in a microlensing event 
by measuring the total equivalent widths of spectral lines on medium-resolution 
spectra with high S/N (as shown in Figure 3), or directly — on high-resolution 
spectra. The change of the Ha line seen in Fig. 3, reflects the expected change 
in the line at different viewing angles — strong at disk center (lens position 0.0) 
and weak at disk limb (lens position 1.0). 

The variation of intensity from disk center to disk limb is studied theo- 
retically by computing the contribution function of each transition; a function 
which defines the line forming region in space (i.e. depth in atmosphere) and 
frequency. Here we define (see Magain 1986) the non-LTE contribution function, 
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Figure 2. Line forming regions, represented by the CFs, of the hy- 
drogen Ha and H/3 lines shown in Figs. lab. The contribution functions 
(CFs) are in units of intensity and are angle averaged, a) The hydro- 
gen Ha line is formed higher in the atmosphere overall, especially its 
line core, b) The hydrogen H/3 line is formed deeper and has well de- 
veloped wings. However, due to the background line and continuum 
opacity these line wings will not be visible at small viewing angles. 
Both lines will have strong, but different, center-to-limb variation. 

CF, for the relative line depression, as: 



where tq is the optical depth at a reference wavelength Ao, tr is defined in terms 
of kr = ki + k c S c /I c kq is the absorption coefficient at Ao, hi and k c are the line 
and continuum absorptions, respectively, S; and S c are the line and continuum 
source functions, respectively, I c is the emergent continuous intensity (if the line 
were absent), and fi = cos 9. The traditional dependence on frequency has 
been omitted from our notation for simplicity. The CF for a given synthesized 
spectral line will indicate where the line is formed, i.e. the line forming region, 
and consequently, at what velocity. We show the CFs for two hydrogen lines in 
our M95-30 model in Figure 2. 

Microlensing suffers from some of the same degeneracies inherent to Doppler 
imaging, namely with respect to the lens path. It could be overcome in the 
coming years with microlensing parallax observations (using spacecraft on the 
Earth orbit). However, unlike Doppler imaging, line profile changes are easy 
to measure and interpret. Of course, the rate of rotation is of no importance 
to microlens imaging - in fact, microlensing can be used to determine vsini 
independently (Maoz & Gould 1994) from the Doppler shifts of the spectral 
lines. 

In order to explore the complete optical spectrum of a lensed star with all 
necessary detail, one needs a very fast code which "microlenses" all frequen- 
cies/wavelengths of an extended source as a function of time. Such an efficient 
method was developed by Heyrovsky & Loeb (1997); Heyrovsky, Loeb, & Sas- 
selov (1997), and is used in the calculations of Fig.3. 

By now there are more than 150 microlensing events detected towards the 
Galactic bulge by the MACHO and OGLE2 projects. The standard light curve 
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of an event is characterized by two observables: its peak amplitude and duration 
(Paczynski 1986). Finite-size effects arise from resolved features on the surface 
of the lensed star. Finite-size effects have been studied as methods to partially 
remove the degeneracy of microlensing light curves through the alterations of the 
standard light curve (Nemiroff & Wickramasinghe 1994; Gould 1994, 1995; Witt 
& Mao 1994; Gould & Welch 1996), its polarized emission (Simmons, Newsam, 
& Willis 1995), spectral shifts due to stellar rotation (Maoz &i Gould 1994), and 
narrow-band photometry in resonance lines (Loeb & Sasselov 1995). Here we 
want to put the emphasis on the inverse problem, that of stellar microscopy — 
reconstructing the stellar surface features and probing the stellar atmosphere. 

Red giants comprise a major component of the stars that are being mon- 
itored in the bulge by the microlensing surveys. For example, in the OGLE 
dataset about 19% of all the stars have radii *S 5R Q (Loeb & Sasselov 1995). 
Therefore, for D QL = 6 kpc, the fraction of giant lensing events with b < R s 
(when the lens transits across the disk of the source star) is of order (R s ) ~ 
4% x {M/0.1M Q y^ 2 . This is a large fraction, and observations are addition- 
ally aided by the intrinsic brightness of the giants. 
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Figure 3. The theoretical model variation of the Ha line (solid line) 
compared to the observed variation (equivalent width) in M95-30 by 
Alcock et al. (1997). Unfortunately, there were no observations made 
near the limb crossing nights (lens position 1.0). The error bars may 
be an underestimate, given the quality of the spectra and the con- 
fused spectral reagion around the Ha line. The observations between 
lens positions 1.4 and 4.0 were obtained on 8 consecutive nights (from 
Heyrovsky, Sasselov, & Loeb 1998). 
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Stellar microscopy can constrain dramatically models of giant atmospheres 
by providing center-to-limb variations of different line profiles in the spectra of 
bulge red giants. Gravitational microlensing offers an easily accessible, imme- 
diate, and inexpensive way to image at least some types of stars. It also offers 
access to stellar populations (in the Galactic bulge and Magellanic Clouds), 
which are beyond the reach of any interferometer. Gravitational lensing does 
resolve, as well as amplify, thus making 17 mag bulge giants as bright as 12 
mag — a target easily accessible to echelle spectrographs. 

There is already a proof-of-concept result, based on the spectroscopy ob- 
tained by the MACHO-GMAN collaboration (Alcock et al. 1997) of event M95- 
30. Our 1-D model of the lensed red giant in that event reproduces well the 
observations of the Ha line shown in Figure 3, as well as the TiO bands and the 
photometry (Heyrovsky, Sasselov, & Loeb 1998). We are currently working on 
the high-quality spectra of M95-30 obtained by Alcock et al. (1997) to improve 
the model. 

3.2. Challenges in Modelling of Cool Giants 

Red giants are a major component of any galactic stellar population. The param- 
eters of the red giant branch are derived from multi-color photometry interpreted 
with stellar atmosphere models (e.g. Vandenberg Sz Bell 1985; Kurucz 1992). 
Such model atmospheres are one-dimensional integrations in semi-infinite slabs 
in hydrostatic and radiative equilibrium. Many problems and the need to im- 
prove these models were reviewed recently by Kurucz (1997). Many assumptions 
and simplifications are being made — some of them have not been verified by 
direct observations. Such examples are the treatment of convective transport 
and small-scale velocity fields (e.g. micro- and macro-turbulence). 

In traditional one-dimensional radiative transfer modeling one usually pa- 
rameterizes small-scale dynamic motions with an assumed microturbulent ve- 
locity distribution. Small-scale in this case refers to motions that occur within 
a photon mean free pathlength, so that these motions contribute to the intrin- 
sic broadening of the line profile. On top of these small-scale motions there 
are motions on larger scales that have to be accounted for in observations with 
limited spatial resolution. These are often called macroturbulence, which is a 
mis-characterization, since most of these motions occur in well organized forms 
like radial- and non-radial pulsations, and convective flows. Macroturbulence is 
applied as additional broadening to the emergent line profiles. As such it has no 
consequences for the radiative transport and energy balance in the atmosphere. 
Microturbulent broadening, on the other hand, does affect the flow of radiation 
through the stellar atmosphere because they provide additional line opacity. 

Although one-dimensional models exist for stars with a wide variety of ef- 
fective temperatures and gravities (e.g., the Kurucz 1992 models or the Allard 
& Hauschildt 1995 models for later type stars) little modeling has been done of 
the multi-dimensional structure of stars other than the Sun (Dravins 1990), and 
practically none for red giants. The stellar evolution models on the other hand 
are dependent on a number of parameters, which describe poorly understood 
physics, e.g. the treatment of convection near the surface of the star. Using scal- 
ing arguments, Schwarzschild (1975) suggested that red giants may have only a 
small number of very large convection cells on their surface. High-resolution mi- 
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crolensing data can provide insights into these issues and help develop efficient 
and physically consistent recipes for calculating small-scale velocity fields. For 
example, assuming that small-scale velocity fields are well represented by a tur- 
bulent spectrum one could model the depth-dependent microturbulent velocity 
with a simple self-consistent model of turbulence in the spirit of Canuto, Gold- 
man, & Mazzitelli (1996), and treat the subgrid cascade. It is clear, however, 
that a 3-D transfer calculation is necessary to provide the initial spectral scales 
for such a scheme (Sasselov 1997b). 

Problems such as described above have hindered significant improvement in 
the determination of fundamental stellar properties over the past 20 years. Just 
to illustrate one of the current problems involving red giants calibrations, con- 
sider the "age of the Universe" debate. Calculating isochrones for globular clus- 
ters — using the same evolutionary models and assumptions, but different model 
atmospheres for T e g and bolometric correction calibrations (Bell-Gustafsson vs. 
Kurucz), leads to 20% difference in age estimates (Salaris et al. 1997; Van- 
denberg 1997). Given the conflicting results in interpreting the HIPPARCOS 
satellite data on subdwarfs calibrations of globular cluster sequences (Reid 1997; 
Pont 1997), the issue of isochrone calibration remains open and very important. 
The need for 2-D and 3-D tools for atmosphere modelling comes also from recent 
successful direct stellar imaging with the Hubble Space Telescope (e.g. Uiten- 
broek, Dupree, & Gilliland 1997). 



4. Conclusions 



Gravitational microlensing offers unique ways to: (1) study depth dependence 
in red giants atmospheres of line formation, small-scale velocity fields, and con- 
vective motions; (2) estimate occurence of surface features like spots and active 
regions in red giants; and (3) derive rotation (-usim) and limb darkening in the 
bulge population of red giants. The imaging of spots and active regions suffers 
from degeneracies similar to those in Doppler imaging. However the ability of 
microlens imaging to resolve center-to-limb variations in spectral lines is excep- 
tional and unique. So is the value of such information to the construction of a 
new generation of giants model atmospheres. 

It is a pleasure to acknowledge my collaborators - D. Heyrovsky and A. 
Loeb, and to thank I. Shapiro for his support of this work. 
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